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. Abstract. We study the large-scale structure formation in the Universe in the frame of scalar- 

' tensor theories as an alternative to general relativity. We review briefly the Newtonian limit of 

non-minimally coupled scalar-tensor theories and the evolution equations of the iV-body system 
i-C ■ that is appropriate to study large-scale structure formation in the Universe. We compute the 

power-spectrum of the universe at present epoch and show how the large-scale structure depends 
Q . on the scalar field contribution. 
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1. Introduction 

During the last two decades high precision observations have established the existence of two 
component of the universe. These are known as dark energy and dark matter. What are their 
fundamental nature? this one of the questions we need to answer [T]. One of the most succesful 
models to explain observations is the ACDM but it still has some problems. 
. In the past we have proposed a dark matter model in the framework of the Newtonian limit of 

a scalar-tensor theory [2J. We have done A^-body simulations using a numerical code developed 
to take into account the scalar field (SF) force contribution which is of Yukawa type in its 
Newtonian limit and we have use it to study several situations, like collision of protogalaxiespJ], 
dynamics and collision of galaxies [H [5], and to study cosmological simulations [6j [8]. Recently 
^ ■ we have used it to address a possible solution to two of the main problems the ACDM has, 

i.e., the tendency to over populate the formed halos of galaxies with satellites and the flatness 
problem [TOKn]. 

Several authors have also been using similar dark matter models to study large-scale structure 
formation, see for instance [12 | fl3j PH] . 

One of the tools to characterize the clustering properties of a universe model is the power 
spectrum. Constraints on cosmological models and on cosmological parameters have been 
derived from measuring power spectrum. Therefore, in this work we consider a dark matter 
model build from the Newtonian limit of a scalar-tensor theory [2] and use it to study the large- 
scale structure formation of the universe by analyzing the power-spectrum versus the parameters 
of the model. 



2. Evolution equations for a ACMD universe within the framework of a 
scalar-tensor theory 

Giving that we are interested in the study of the large-scale structure formation of the universe 
we will use the Newtonian limit of the scalar-tensor theory This limit applies in a region small 
compared to the Hubble length cH^ 1 , with c the speed of light and H the hubble factor, and 
large compared to the Schwarzschild radii of any collapsed objects. 

In this case the Vlasov-Poisson equation in an expanding universe describes the evolution of 
the six-dimensional, one-particle distribution function, /(x, p)[15j. The Vlasov equation is, 

— + j.- mV$jv(x) ■ — = (1 

at ma z ax op 

where x is the comoving coordinate, p = mo 2 x, m is the particle mass, and <£at is the self- 
consistent gravitational potential given by the Poisson equation, 

V 2 $iv(x) = 4ttG n a 2 [p(x) - p b (t)}} (2) 

where pb is the background mass density. Eqs. ([1]) and ([2]) form the Vlasov-Poisson equation, 
constitutes a collisionless, mean-field approximation to the evolution of the full iV-body 
distribution. An A^-body code attempts to solve Eqs. ([TJ and ([2]) by representing the one- 
particle distribution function as 

N 

/(x,p)=^<5(x-x i )5(p-p i ) (3) 

i=i 

Substitution of (|3|) in the Vlasov-Poisson system of equations yields the exact Newton's equations 
for a system of N gravitating particles. 

In the Newtonian limit of STT of gravity, the Newtonian motion equation for a particle i is 
written as [7] 

x, + 2 x, = E T^"^ ^(1* " x ^'l' «> A ) ( 4 ) 

a 6 1 + a rr! x» — x,- r 

where the sum includes all periodic images of particle j, and FsF(r, a, A) is 

F SF (r,a,X) =l + a (l + £) ^ r/A (5) 

which, for small distances compared to A, is FsF(r <C A, a, A) « 1 + a (1 + v) and, for long 
distances, is Fsf(t 3> A, a, A) w 1, as in Newtonian physics. 

To simulate cosmological systems, the expansion of the universe has to be taken into account. 
Also, to determine the nature of the cosmological model we need to determine the composition 
of the universe, i. e., we need to give the values of f2j for each component i, taking into account 
in this way all forms of energy densities that exist at present. 

In this work we will consider a model with only two energy density contribution. One which is 
a pressureless and nonbaryonic dark matter with £Idm ~ 0.3 that does not couple with radiation. 
Other, that will be a cosmological constant contribution ~ 0.7. Also, the model must be 
consistent with a static SF according to the Newtonian limit we are using in this work. Thus, 
the scale factor, a(t), is given by the following Friedman model [8 , 

fiflMo + ^ao a 3 . f, ^DMO + ^A0 \ 
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where H = a/ a, Qdmo and Oao are the matter and energy density evaluated at present, 
respectively. We have also demanded that the universe be a flat one. To be consistent with 
the CMB spectrum and structure formation numerical experiments, cosmological constraints 
must be applied on a in order for it to be within the range (—1,1) [161 H3 fTHl f!9] . 
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Figure 1. Power spectrum of our scalar field universe model without and with SF. See text for 
details. 

3. Results 

In this section, we present results of cosmological simulations of a ACDM universe with and 
without SF contribution. We start our simulations with an initial distribution of N = 2 x 32 3 
particles in a box with sides of L = 50/t -1 Mpc at z = 10. This case is similar to the one that 
comes with Gadget 1 [20]. At present epoch, Q m Q = 0.3, f^o = 0.7, Hq = WOh km/s/Mpc, 
h = 0.7. We restrict the values of a to the interval (-1, 1) P3 El US US] and use A = 5 Mpc, 
since this scale turns out to be an intermediate scale between the size of the clump groups and 
the separation of the formed groups. However, we have analyzed elsewhere several values of 
A[8]. Other universe models have been analyzed in [S1E1IS]- 

Because the visible component is the smaller one and given our interest to test the 
consequences of including a SF contribution to the evolution equations, our model excludes 
gas particles, but all its mass has been added to the dark matter. In Fig. 1 we show the power 
spectrum at redshift z = of our ACDM model for several values of a. The line without symbols 
shows the case of a ACDM model without the contribution of the scalar field. The line with 
(+) is the case of a ACDM model with scalar field with a = 1, the line with (*) uses a = —1/2, 
and the line with (o) uses a = —1/4. One notes clearly how the SF modifies the clustering of 
matter of the system. The most dramatic cases are the models with a = —1/2 and a = 1. 

We now analyze the general effect the parameter a has on the clustering. For a given A 
the role of a in our approach can be seen in equations (H])-© and is contained in the factor 
Fsf(i", a, A). The factor F$f augments (diminishes) for positive (negative) values of a for small 
distances compared to A, resulting in more (less) structure formation for positive (negative) 
values of a compared to the ACDM model. In the case of the line with (+), for r <C A, the 
effective gravitational pull has been augmented by a factor of 2, in contrast to case shown with 
line and (*) where it has diminished by a factor of 1/2; in the model with the line and (o) the 
pull diminishes only by a factor of 3/4. That is why one observes for r < A more structure 
formation in case with line and (+), less in case with line and (*), and lesser in model with line 
and (o). The effect is then, for a growing positive a, to speed up the growth of perturbations, 
then of halos and then of clusters, whereas negative a values (a — > —1) tend to slow down the 
growth. 



4. Conclusions 

The general gravitational effect of the scalar-tensor dark matter model is that the interaction 
changes by a factor Fsf(t, ol, A) ~ 1 + cu(l + r/A) for r < A in comparison with the Newtonian 
case. For negative values of a the effect is to diminish the formation of structures and the 
opposite ocurrs for positive values of a. For r > A the dynamics is essentially Newtonian. 
By computing the power spectrum we were able to determine the effects of the scalar field 
contribution to the galaxy clustering. Of course we will need to do a systematic and with 
greater resolution study and make detailed comparisons with observations to support or to rule 
out a dark model based on scalar-tensor theories. 
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